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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SJN 
UNITS  OF  MEASUREMENT 

U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con¬ 
verted  to  metric  (SI)  units  as  follows: 


_ Multiply _ 

acres 

acre-feet 

Btu  (International  Table) 

Btu  (International  Table)/ 
lb/°F 

cubic  feet  per  second 
Fahrenheit  degrees 
feet 

feet  per  second  per  second 


_ lz _ 

4046. 856 
1233.482 
1055.056 
4186.800 

0.02831685 

5/9 

0.3048 

0.3048 


miles  (U.  S.  statute)  1.609344 

pounds  (mass)  per  cubic  foot  16.01846 
square  feet  0.09290304 

tons  (2,000  lb,  mass)  907.1847 


_ To  Obtain _ 

square  metres 
cubic  metres 
joules 

joules  per  kilogram  kelvin 

cubic  metres  per  second 
Celsius  degrees  ox  Kelvins* 
metres 

metres  per  second  per 
second 

kilometres 

kilograms  per  cubic  metre 
square  metres 
kilograms 


*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (i?)  read¬ 
ings,  use  the  following  formula:  C  *  (5/9) (F  -  32).  To  obtain  Kelvin 
(K)  readings,  use:  K  -  (5/9) (F  -  32)  +  273.15. 
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EVALUATION  OF  IN- RESERVOIR  COFFERDAM  ON  RICHARD  B.  RUSSELL 
RESERVOIR  AND  HYDROPOWER  RELEASES 

Hybrid  Model  Invest igat ion 

PART  Is  BACKGROUND 

1.  A  coldwater  release  (less  than  70°K)*  with  a  dissolved  oxygen 
(D.O. )  concentration  of  6  mg/2.  is  desired  from  the  proposed  Richard  B. 
Russell  Ft- servo ir  (RBRR)  and  pumped  storage  hydropower  facilities.  Dur¬ 
ing  stratification  of  the  reservoir,  the  respective  vertical  distribu¬ 
tions  of  temperature  and  D.O.  essentially  preclude  withdrawal  and  release 
of  water  that  meets  both  objectives.  If  sufficient  epilimnetlc  water 
is  withdrawn  to  meet  the  D.O.  objective,  the  upper  bound  of  the  tempera¬ 
ture  criterion  is  exceeded.  Conversely,  if  adequate  cold  hypclimnetic 
water  is  withdrawn  to  meet  the  temperature  objective,  the  withdrawn  D.O. 
concentration  is  less  than  that  desired j  thus,  if  adequate  reaeration 
does  not  occur  in  or  below  the  hydropower  structure,  the  D.O.  objective 
will  not  be  met.  To  circumvent  this  potential  problem,  the  U.  S.  Army 
Engineer  District,  Savannah  (SAS)  will  (a)  place  the  center  line  of  the 
hydropower  intakes  approximately  70  ft  below  the  surface  of  the  pool  to 
withdraw  cold  water  and  (b)  install  an  innovative  oxygenation  system  in 
RBRR  to  increase  the  D.O.  concentration  in  the  withdrawal  zone  of  the 
reservoir.  Obviously,  large  quantities  of  molecular  oxygen  are  involved 
in  this  procedure.  As  a  result,  SAS  and  the  U.  S.  Army  Engineer  Divi¬ 
sion,  South  Atlantic  (SAD),  requested  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  to  investigate  the  relative  effectiveness 
and/or  feasibility  of  various  in-reservoir  structural  alternatives  in 
enhancing  the  in  situ  D.O.  in  tiie  reservoir  withdrawal  zone,  thus  reduc¬ 
ing  the  quantity  (or  rate)  of  oxygen  required  for  injection.  The  in¬ 
reservoir  design  options  identified  by  SAS  and  SAD  were: 

*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  units  is  presented  on  page  3. 


a.  Option  A;  removal  of  the  cofferdam  (as  originally 
planned) . 

b.  Option  B:  continuous  cofferdam  as-constructed  with  a 
crest  at  el  378.*  Diversion  breach  would  be  filled. 

c. .  Option  C:  cofferdam  left  in  place  as-constructed  and 

without  filling  the  second  stage  50-ft-long  diversion 
breach  (bottom  at  el  310). 

d.  Option  D:  a  sheet-pile  cutoff  wall  installed  on  the  crest 
of  the  continuous  cofferdam  with  the  top  at  el  400, 

A  schematic  of  the  geometries  is  presented  in  Figure  1. 

Project  Description 

2.  The  RBRR  is  a  multipurpose  project  currently  nearing  comple¬ 
tion  of  construction  by  SAS.  The  ^reject,  which  will  impound  the 
Savannah  River  between  the  headwaters  of  Clarks  Hill  Reservoir  and  Hart¬ 
well  Dam  (Figure  2),  is  designed  to  provide  hydropower.  The  dam  con¬ 
sists  of  a  195-ft-high,  1884-f t-long  concrete  gravity  structure  flanked 
by  two  earth  embankments.  The  powerhouse  will  house  four  75-MW  conven¬ 
tional  units  and  four  75-MW  pump  units.  Total  capacity  of  the  system 

is  600  MW.  Withdrawal  for  these  units  will  be  accomplished  through 
eight  penstocks;  each  penstock  has  a  cross-section  area  of  approximately 
5400  ft  and  is  located  at  center-line  el  393. 

3.  The  RBRR  will  have  a  surface  area  of  26,650  acres  when  filled 
to  maximum  power  pool,  el  475.  At  this  elevation  the  reservoir  will 
have  a  storage  volume  of  approximately  1,026,200  acre-ft  and  will  extend 
up  the  Savannah  River  to  the  vicinity  of  the  Hartwell  Dam.  Depth  of  the 
reservoir  at  maximum  power  pool  will  be  165.0  ft.  The  release  of  exces¬ 
sive  flood  flows  will  be  accomplished  by  a  500-ft-long  ogee  spillway. 

Approach 

4.  The  hybrid  modeling  approach  was  used  in  this  study.  A 


*  All  elevations  (el)  cited  herein  are  in  feet  referred  to  mean  sea 
level  (rasl) . 


ft  475 


'RULE  CURVE 


cl 


ft  m 


ILSl 


CL  310 


a.  Cofferdam  removed,  Option  A 


■rule  cur  vc 


EL  393  _ 


b.  Continuous  or  breached  cofferdam  (Options  B  and  C) 


» RULE  CURVE 


EL  437 


EL  400 


EL  393 _ _ 

CL  370 


1100  FT 


c.  Raised  cofferdam.  Option  D 

Figure  1.  Schematic  of  in-reservoir  geometrical  alternatives 


. .  ••  •  •-  - 


Figure  2.  Location  map 

numerical  simulation  model,  WESTEX,  was  adopted  from  previous  analyses 
of  pumped  storage  projects,*  Details  of  the  numerical  model  are  pro¬ 
vided  in  Appendix  A,  Algorithms  that  addressed  the  selective-withdrawal 
and  purapback  processes  were  developed  and  calibrated  with  data  obtained 
from  an  undlstorted-scale  physical  model  that  simulated  the  hydrodynamics 
occurring  near  the  dam.  The  physical  model  was  constructed  to  a  1:150 
scale  and  reproduced  the  dam  and  approximately  3000  ft  of  the  upstream 


*  M.  S.  Dortch  el  al.  1976  (Dec).  'Dickey-Lincoln  School  Lakes 
Hydrothermal  Model  Study;  Hydraulic  Model  Investigation,"  Technical 
Report  H-76-22,  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE, 
Vicksburg,  Miss. 

D.  G.  Fontane  et  al.  1977  (Apr).  "Marysville  Lake  Hydrothermal 
Study:  900-MW  Project;  Hydraulic  and  Mathematical  Model  Investiga¬ 
tion,"  Technical  Report  H-77-5,  Report  1,  U.  S.  Army  Engineer  Water¬ 
ways  Experiment  Station,  CE,  Vicksburg,  Miss. 
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topography.  Additional  details  of  the  physical  model  and  the  selective 
withdrawal  and  the  pumpback  descriptions  developed  from  the  model  data 
are  provided  in  Appendix  B.  The  numerical  model  was  used  with  histori¬ 
cal  hydrology  and  meteorology  for  1955,  1958,  1966,  and  1967  to  simulate 
the  reservoir  temperature  and  vertical  distributions  of  D.O.  and  the 
release  temperature  and  D.O.  concentrations  expected  during  the  four 
study  years  with  each  of  the  in-reservoir  design  alternatives.  The  re¬ 
spective  predictions  for  each  in-reservoir  design  alternative  were  then 
compared  to  ascertain  their  relative  merits  in  enhancing  the  D.O.  dis¬ 
tribution  in  the  reservoir  and  the  D.O.  concentration  in  the  release. 
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PART  II:  ACQUISITION  AND  DEVELOPMENT  OF  DATA 


5.  The  simulation  model  (WESTEX)  required  data  for  lake  inflows, 
outflows,  inflow  stream  temperature,  various  meteorological  data,  and 

a  schedule  reflecting  daily  power  (generation  and  pumpback)  operations. 
The  WESTEX  model  also  required  determination  of  coefficients  for  surface 
heat  exchange  and  for  internal  mixing.  Various  data  required  by  the 
simulation  model  are  described  below. 

Meteorology 

6.  Meteorological  data  for  the  Atlanta,  Georgia,  Class  A  weather 
station  was  obtained  from  the  National  Climatic  Center,  Asheville,  North 
Carolina.  The  required  data  consisted  of  dry  bulb  and  dewpoint  tempera¬ 
tures,  wind  speed,  and  cloud  cover.  From  these  data,  daily  average 
values  of  equilibrium  temperature,  surface  heat  exchange  coefficient, 
and  net  solar  radiation  quantity  were  computed.  For  each  study  year, 
these  parameters  were  computed  by  the  U.  S.  Army  Engineer  Division, 

Ohio  River. 


Lake  Geometry 

7.  The  elevation-capacity  data  for  RBRR  and  the  location  and 
design  of  the  penstocks  were  supplied  by  SAS. 

Inflow  Quantity  and  Temperature 

8.  Both  the  inflow  quantities  and  temperature  were  obtained  from 
observed  outflow  quantities  and  temperatures  from  Hartwell  Dam,  which 
will  release  directly  into  RBRR.  These  values  were  supplied  by  SAS. 

Power  Operations 

9.  Conventional  power  generation  releases  at  RBRR  will  range  from 
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0  to  35,000  cfs  during  a  one-day  period.  Pumped  storage  operations  at 
RBRR  will  increase  the  upper  bound  of  the  hydropower  releases  to 
60,000  cfs.  The  proposed  operation  schedule  '.ndicates  that  the  pumpback 
rates  from  the  headwaters  of  Clarks  Hill  to  RBRR  will  remain  nearly  con¬ 
stant  at  20,000  cfs.  The  load  curve  used  by  Fontane  and  Bohan*  appears 
in  Figure  3  and  was  used  for  this  study  as  directed  by  SAS.  For  each 
simulation  day,  a  schedule  of  operations  representing  one  or  two  genera¬ 
tion  cycles  and  one  (or  none,  if  conventional  power  only)  pumpback  cycle 
was  modeled  after  Figure  3.  The  number  of  hours  of  generation  was  held 
constant  for  each  day  as  shown  in  the  load  curve.  The  generation 
release  rates  matched  the  outflow  data  supplied  by  SAS. 


2  *o, 
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Figure  3.  RBRR  load  curve 


Coefficient  Selection 


10.  As  has  been  discussed  previously,  the  WESTEX  model  requires 
determination  of  coefficients  of  surface  heat  exchange  and  internal 
mixing.  Two  surface  heat  exchange  coefficients  must  be  selected:  B  , 
the  shortwave  radiation  absorbed  in  the  surface  layer,  and  X  ,  the 
light  extinction  coefficient.  Similarly,  two  mixing  coefficients, 
and  ,  must  be  selected.  These  coefficients  are  analogous  to  mixing 
produced  by  vertical  diffusion.  The  four  coefficients  were  predicated 
upon  a  comparison  of  predicted  profiles  for  RBRR  without  pumpback  with 
observed  profiles  from  Hartwell  Reservoir.  Hartwell  Reservoir  is 


*  D.  G.  Fontane  and  J.  P.  Bohan.  1974  (Dec).  "Richard  B.  Russell  Lake 
Water  Quality  Investigation;  Hydraulic  laboratory  Investigation,” 
Technical  Report  H-74-14,  U.  S.  Army  Engineer  Waterways  Experiment 
Station,  CE,  Vicksburg,  Miss. 


10 


PART  III:  RESULTS 


11.  The  relative  impact  of  the  respective  geometrical  configura¬ 
tions  upon  the  in-reservoir  and  release  water  quality  can  be  demon¬ 
strated  by  comparing  representative  in  situ  and  release  water  quality 
predictions.  The  in-reservoir  results  trill  be  discussed  first. 

In-Reservoir  Dissolved  Oxygen  and  Temperature 

Dissolved  oxygen 

12.  Predicted  D.O.  profiles  for  the  respective  in-reservoir  geome¬ 
tries  are  presented  in  Figures  4  and  5.  Although  these  profiles  are  for 
only  27  September  of  the  various  study  years,  they  are  indicative  of 

the  relative  effects  that  each  geometry  has  upon  the  D.O.  distribution 
during  the  late  phases  of  stratification.  As  Indicated  in  Figures  4 
and  5,  the  most  favorable  D.O.  distribution  during  the  late  phases  of 
stratification  is  obtained  if  the  cofferdam  is  removed  (Option  A, 

Figure  1).  With  the  cofferdam  removed,  a  larger  volume  of  water  is 
entrained  from  the  hypollmnlon  into  the  pumpback  jet.  After  entrainment, 
the  low  D.O.  hypollnmetic  water  mixes  with  metalimnetic  (or  eplllmnetlc) 
waters.  Simultaneously,  water  with  higher  D.O.  from  the  epllimnion  (or 
metalimnion)  moves  downward.  Since  oxygen  is  continually  supplied  to 
the  epllimnion  through  the  air-water  interface,  and  since  the  volume 
from  the  hypolinmlon  pumped  during  a  pumpback  operation  is  small  com¬ 
pared  with  the  total  volume  of  the  epllimnion  or  metalimnion,  signifi¬ 
cant  D.O.  depletion  in  the  upper  region  does  not  occur.  The  continual 
supply  of  high  D.O.  water  and  the  redistribution  of  this  water  results  in 
a  tendency  to  improve  the  D.O.  distribution  in  the  reservoir. 

13.  Any  horizontal  obstruction  in  the  reservoir  that  reduces  the 
quantity  of  water  removed  from  the  hypolimnion  by  entrainment  will  tend 
to  decrease  vertical  mixing  during  pumpback.  With  the  continuous  and 
breached  cofferdams  (Options  B  and  C,  Figure  1),  the  physical  model 
indicated  13  percent  less  entrainment  from  the  hypolimnion  than  would 
occur  if  the  cofferdam  were  removed.  Thus,  a  less  favorable  D.O. 
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Figure  5.  Predicted  D.O.  distribution  with  the  cofferdam 
removed  and  with  the  raised  cofferdam 

indicated  transient  flow  through  the  breach  between  pumpback  and  with¬ 
drawal  phases.  As  a  result,  if  the  cofferdam  is  breached  and  poor  water 
quality  exists  in  the  lower  layers,  gravity  flow  through  the  breach 
during  the  time  interval  between  pumpback  and  release  could  result  in 
poor  water  quality  between  the  dam  and  the  cofferdam.  During  subsequent 
hydropower  release,  a  pulse  of  low  D.O.  water  could  occur, 

15.  Placing  a  sheet-pile  cutoff  wall  with  a  crest  at  el  400  (Op¬ 
tion  D,  Figure  1)  severely  mitigates  vertical  mixing  during  pumpback. 
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As  compared  with  removal  of  the  cofferdam  (Option  A),  90  percent  leee 
entrainment  from  the  hypollmnion  was  measured  In  the  physical  model  with 
the  raised  cofferdam.  Depending  upon  the  simulation  year,  anaerobic 
conditions  developed  in  the  hypollmnlon  in  late  June  to  mid-July.  By 
the  end  of  September,  the  thickness  of  the  anaerobic  sone  had  Increased 
to  between  65  and  80  ft  depending  upon  the  simulation  year.  As  a  result 
this  alternative  is  not  considered  viable  and  will  not  receive  further 
consideration. 

16.  A  synopsis  of  the  in-reservoir  analysis  with  respect  to  the 
D.O.  distribution  is  presented  in  Table  1. 

In-reservoir  temperature 

17.  Representative  temperature  profiles  are  presented  in  Figure  6 
Temperature  distributions  produced  in  the  reservoir  by  withdrawal  over 
the  cofferdam  were  colder  than  those  predicted  with  the  cofferdam  re¬ 
moved.  The  withdrawal  distribution  is  higher  in  the  reservoir  as  a 
result  of  withdrawal  over  the  cofferdam  which  results  in  the  release  of 
warmer  water  and  the  retention  of  colder  water. 

Release  Dissolved  Oxygen 

18.  The  predicted  release  D.O.  concentration  as  a  function  of 
time  for  the  1966  and  1967  study  years  is  presented  in  Figure  7.  These 
results  are  similar  to  the  results  obtained  for  the  1955  and  1958  study 
years  and  may  be  used  to  demonstrate  general  findings. 

Qualitative  observations 

19.  Prior  to  June  and  after  September  (or  October  for  conditions 
similar  to  1955),  the  predicted  release  D.O.  concentrations  with  the 
respective  cofferdam  geometries  were  essentially  equivalent.  The  differ 
ence  in  the  release  D.O.  for  the  respective  geometries  during  these 
periods  was  less  than  0.25  mg/ JU  Thus,  if  the  reservoir  was  unstrati¬ 
fied  or  only  weakly  stratified,  the  in-reservoir  geometry  did  not  sig¬ 
nificantly  impact  the  release  D.O. 

20.  The  release  D.O.  concentration  was  more  sensitive  to  the  ln- 
reservoir  geometry  during  significant  stratification.  In  general,  the 


15 


HEIGHT  FROM  BOTTOM.  FT 


200 

160 

120 

80 

40 

0 

200 

160 

120 

80 

40 

0 


■ 

■ 

■ 

S3 

■ 

■ 

■ 

■ 

*2 

m 

■ 

■ 

■ 

BS 

■ 

m 

■ 

■ 

n 

■ 

■ 

■ 

■ 

ii 

■ 

■ 

■ 

TEMPERATURE,  °C  TEMPERATURE,  °C 


LEGEND 

- COFFERDAM  REMOVED 

CONTINUOUS  OR  BREACHED 
COFFERDAM 


Figure  6.  Predicted  temperature  distribution  without  the  cofferdam 
and  with  the  continuous  or  breached  cofferdam 

release  D.O.  concentration  with  the  cofferdam  either  continuous  or 
breached  tended  to  be  superior  to  those  predicted  with  the  cofferdam 
removed.  These  qualitative  results  are  indicated  in  Figure  7.  Essen¬ 
tially  identical  predictions  were  obtained  for  the  breached  (Option  C) 
and  continuous  cofferdam  (Option  B)  for  the  reasons  indicated  in  para¬ 
graph  14.  These  qualitative  results  reflect  the  general  trends  but 
inherently  do  not  quantify  the  relative  merits  of  the  respective  geomet¬ 
rical  configuration. 
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Quantitative  results 

21.  Deviation  plots  for  each  simulation  year  are  presented  in 
Figure  8.  The  ordinate  is  the  predicted  release  D.O.  concentration 
with  the  cofferdam  (Option  B  or  C)  minus  the  predicted  D.O.  with  the 
cofferdam  removed  (Option  A).  A  positive  deviation  indicates  that  the 
release  D.O.  concentration  with  the  cofferdam  would  be  greater  by  the 
indicated  magnitude  than  that  concentration  release  without  the 
cofferdam. 

22.  Between  1  June  and  mid-July ,  the  cofferdam  resulted  in  slight 
improvement  in  the  release  L.O.  as  compered  with  that  with  the  cofferdam 
removed.  The  magnitude  of  enhancement  was  somewhat  erratic  as  a  result 
of  the  dependency  upon  the  hydrological  and  meteorological  conditions; 
however,  in  all  cases  maximum  Improvement  was  less  than  1  mg/ 2  and  the 
average  improvement  over  this  time  frame  was  no  more  than  0.7  mg/2.. 
Beginning  in  late  July  to  mid-August,  a  period  existed  for  each  study 
year  in  which  the  D.O.  enhancement  was  in  excess  of  1  mg/ 2.  This  period 
varied  from  25  to  58  days  (depcndin3  on  the  study  year)  and  the  magni¬ 
tude  of  D.O.  enhancement  generally  varied  between  1  mg/2  and  1.8  mg/2. 
Onl_’  during  1966  did  D.O.  enhancement  exceed  2  rng/f  and  this  was  for 
only  a  three-day  perio  t.  A«?  Indicated  in  Figure  8,  the  period  of  maxi¬ 
mum  D.O.  enhancement  dominates  the  deviation  plot  and  the  magnitude  of 
D.O.  enhancement  usually  drops  rapidly  on  both  sides  of  this  region. 

Idealized  predictions  of 
required  injected  oxygen 

23.  To  wichdrai  water  from  the  reservoir  that  has  an  average  D.O. 
concentration  of  6  mg/ 2,  oxygen  would  have  to  be  added  to  the  withdrawal 
zone  in  the  reservoir  during  periods  in  which  the  volumetrically 
weighted  average  D.O.  concentration  in  the  release  is  less  than  6  mg/ 2. 
Typically,  the  release  D.O.  concentration  was  less  than  6  mg/2  between 
mid-May  and  October  (Figure  7).  The  idealized  daily  average  injection 
rate  required  to  eradicate  deficits  below  6  mg/t  was  computed  by  multi¬ 
plying  the  average  daily  release  by  the  volumetrically  weighted  average 
D.O.  below  6  mg/2.  Predicted  daily  average  injection  rates  for  1966  are 
presented  in  Figure  9.  Although  the  predicted  injection  histories  are 
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unique  for  each  study  year,  these  resultB  are  representative  of  the 
results  for  all  simulation  years. 

24.  The  temporal  injection  rate  required  over  a  simulation  de¬ 
pends  upon  a  large  number  of  variables.  However,  at  any  particular 
time  it  is  essentially  dependent  upon  (a)  the  vertical  temperature  and 
D.O.  distributions  and  (b)  the  operational  methodology.  Both  of  these 
effects  are  apparent  in  Figure  9.  During  maximum  stratification,  the 
least  favorable  D.O.  distribution  develops  in  the  metalimnlon  and 
hypolimnion.  The  center  line  of  the  penstock  is  approximately  70  ft 
below  the  surface  and  thus  significant  quantities  of  water  are  withdrawn 
from  the  metalimnlon  and  hypolimnion.  If  large  releases  occur  during 
periods  of  low  D.O. ,  large  injection  rates  are  required.  Conversely,  if 
small  releases  are  made,  smaller  injection  rates  are  required.  Thus, 
the  predicted  injection  rates  reflect  the  coupling  between  operations 
and  the  D.O.  in  the  withdrawal  zone.  Further,  D.O.  concentrations  in 
the  withdrawal  zone  varies  rather  slowly  whereas  the  hydropower  opera¬ 
tions  are  highly  transient.  As  a  result,  erratic  daily  injection  rates 
are  computed. 

25.  The  maximum  daily  injection  rate  (neglecting  spikes  created 
by  the  operational  strategy)  computed  for  each  study  year  is  presented 
in  the  following  tabulation: 

Computed  Maximum  Injection  Rate 


tons  per  day 

Year 

Cofferdam 

Removed 

Continuous 

Cofferdam 

Breached 

Cofferdam 

1955 

110 

62 

62 

1958 

97 

55 

55 

1966 

110 

75 

75 

1967 

95 

70 

70 

26.  A  more  comprehensive  comparison  of  the  respective  injection 
rates  can  be  made  from  the  deviation  plots  for  each  simulation  ye.  r  (Fig¬ 
ure  10),  The  ordinate  is  the  predicted  injection  rate  with  the  ct ffer- 
dam  minus  the  predicted  rate  with  the  cofferdam  removed.  A  negative 
magnitude  Indicates  that  the  required  injection  rate  with  the  cofferdam 
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•  10-  Computed  differences  in  idealized  injections  rates.  Rate* 
cofferdam  (continuous  or  breached)  minus  rate  without  cofferdam 


is  the  indicated  magnitude  less  than  that  (predicted)  without  the 
cofferdam.  As  compared  with  removal,  retention  of  the  cofferdam  con¬ 
sistently  resulted  in  smaller  required  oxygen  injection  rates.  In  gen¬ 
eral,  prior  to  mid-July  and  after  September  the  difference  in  the  pre¬ 
dicted  injection  rate  was  relatively  insignificant.  During  August  and 
September  the  maximum  difference  in  the  predicted  injection  rates  was 
between  35  and  50  tons  per  day  (depending  upon  the  simulation  year). 

27.  The  total  quantity  of  oxygen  required  in  the  prototype  for 
either  (or  both)  in-reservoir  geometries  may  be  significantly  different 
from  those  predicted  by  the  above  idealized  computations.  The  required 
injected  quantity  will  depend  upon  the  efficiency  of  the  particular  in¬ 
jection  system  in  enhancing  the  D.O.  distribution  within  the  unique 
withdrawal  zone  for  the  in-reservoir  geometry  selected.  It  is  improb¬ 
able  that  a  particular  injection  system  and  operational  scheme  will  be 
equally  effective  for  all  in-reservoir  geometries  investigated  since  the 
withdrawal  and  pumpback  velocity  distributions  in  the  reservoir  are 
different  for  various  geometries.  Thus,  from  idealized  computations  it 
is  not  possible  to  accurately  deduce  either  quantitatively  (or  rela¬ 
tively)  the  injected  oxygen  required  in  the  respective  prototype  config¬ 
urations.  These  calculations  serve  as  indicators  based  upon  the 
integrated  deficit  and  an  idealized  injection  system. 

Release  Temperature 

28.  A  representative  temperature  deviation  plot  is  presented  in 
Figure  11.  A  positive  deviation  indicates  that  the  release  temperature 
with  the  cofferdam  would  be  larger  by  the  indicated  magnitude  than  that 
release  without  the  cofferdam.  The  cofferdam  resulted  in  slightly 
warmer  releases  because  of  the  higher  withdrawal  distribution.  The 
temperature  increase  during  the  stratification  season  was  typically  be¬ 
tween  0.5  and  2°C.  The  duration  of  the  increases  was  approximately  six 
to  eight  months  (depending  upon  the  simulation  year). 
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PART  IV:  CONCLUSIONS 


29.  The  relative  effectiveness  and/or  feasibility  of  four  in¬ 
reservoir  structural  alternatives  in  enhancing  reservoir  and  release 
water  quality  was  predicted.  Predicated  upon  analysis  of  these  predic¬ 
tions,  the  alternatives  should  be  ranked  in  the  following  order: 

Rank  _ Option _ 

1  Removal  of  the  cofferdam 

2  Continuous  cofferdam 

3  Breached  cofferdam 

A  Raised  cofferdam 

30.  The  least  favorable  option  is  the  raised  cofferdam.  Even 
though  release  D.O.  concentrations  are  increased  by  raising  the  crest  to 
el  400,  substantially  reduced  water  quality  may  result  in  the  reservoir. 
Reasonable  depletion  rates  would  result  in  a  rather  large  anaerobic  re¬ 
gion  in  the  hypolimnion.  This  occurs  because  the  high  crest  elevation 
minimizes  advection  and  vertical  mixing  at  the  deeper  elevations.  This 
option  is  not  recommended. 

31.  Although  the  predicted  water  quality  (both  in  the  reservoir 
and  in  the  release)  was  basically  the  same  for  the  continuous  and 
breached  cofferdam,  the  latter  option  is  not  considered  desirable.  The 
breach  provides  a  flow  path  below  the  crest  of  the  cofferdam.  Between 
pumpback  and  generation,  gravity  flow  through  the  breach  could  result 
in  lower  quality  water  between  the  cofferdam  and  dam.  During  the  sub¬ 
sequent  hydropower  release,  a  pulse  of  water  with  relatively  low  D.O. 
concentrations  could  be  released.  This  flow  process  was  observed  in  the 
physical  model  but  significant  investigation  of  this  flow  was  outside 
the  scope  of  the  numerical  study. 

32.  The  effect  of  the  cofferdam  (continuous  or  breached)  on  the 
performance  of  the  proposed  RBRR  oxygen  injection  system  is  unknown  be¬ 
cause  the  injection  system  was  developed  from  pilot  studies  at  Clarks 
Hill  Reservoir  which  did  not  have  a  submerged  cofferdam.  Results 

from  the  physical  model  indicated  that  the  cofferdam  significantly 
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affects  the  hydrodynamics  of  both  selective  withdrawal  and  pumpback. 

As  a  result,  the  presence  of  a  cofferdam  could  have  an  adverse  effect 
upon  the  projected  efficiency  of  the  injection  system  in  enhancment  of 
the  D.O.  concentration  in  the  withdrawal  zone. 

33.  Any  horizontal  obstruction  in  the  reservoir  which  reduces  the 
water  withdrawn  from  the  hypolimnion  by  withdrawal  or  pumpback  (entrain¬ 
ment)  has  the  effect  of  decreasing  vertical  mixing  and  advection.  The 
decreased  mixing  and  increased  residence  time  can  result  in  a  degrada¬ 
tion  in  the  D.O.  distribution  (with  time)  in  the  lower  layers  if  oxygen 
demands  exists.  The  impact  of  a  depletion  rate  of  0.1  mgM/day  upon  the 
D.O.  distribution  during  September  of  each  study  year  is  demonstrated  in 
Figure  4.  If  severe  loading  exists  in  the  hypolimnion,  the  impact  could 
be  more  severe. 

34.  Removal  of  the  cofferdam  should  complement  the  oxygen  system. 
The  proposed  oxygen  injection  system  is  designed  to  enhance  the  D.O.  dis¬ 
tribution  in  the  zone  of  withdrawal  such  that  the  flow-weighted  average 
D.O.  concentration  is  equal  to  or  in  excess  of  6  mg/i.  At  any  particu¬ 
lar  time,  the  desired  Increase  in  the  D.O.  distribution  depends  upon  the 
deficit  and  velocity  distribution  within  the  withdrawal  zone.  It  is  dif¬ 
ficult  to  design  an  optimum  variable  rate  system  that  matches  this  ob¬ 
jective  at  all  times  because  the  deficit  distribution  varies  with  time. 
Thus,  the  system  has  to  be  designed  for  the  maximum  anticipated  deficit. 
However,  it  must  also  be  capable  of  distributing  the  oxygen  appropriately 
within  the  zone  of  withdrawal.  The  characteristics  of  the  withdrawal 
zone  (location  of  the  limits  and  velocity  distribution)  vary  with  the 
hydrological  and  meteorological  conditions  and  the  cofferdam  configura¬ 
tion.  Removal  of  the  cofferdam  should  make  designing  an  injection  sys¬ 
tem  less  difficult.  The  enhanced  vertical  mixing  that  results  from  the 
pumpback  process  will  tend  to  distribute  the  D.O.  and  make  placement  of 
oxygen  within  the  withdrawal  zone  less  critical. 

35.  The  deeper  withdrawal  distribution  that  would  result  with  the 
cofferdam  removed  would  provide  slightly  cooler  releases  (1°C  to  2°C) 
than  with  the  other  in-reservoir  structural  alternatives.  During  periods 
of  relatively  strong  stratification,  the  release  temperatures  would  be 
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more  favorable.  If  the  temperature  criterion  is  exceeded,  this  geo¬ 
metrical  alternative  would  tend  to  reduce  its  duration  and  intensity. 
This  should  be  beneficial  in  meeting  the  coldwater  temperature 
objective. 

36.  The  sizing  of  the  oxygen  injection  plant  was  outside  the 
scope  of  the  WES  effort.  Detailed  development  of  the  oxygen  injection 
system  by  SAS  is  described  in  Design  Memorandum  3.  Idealized  injection 
rates  were  included  herein  only  to  serve  as  indicators. 
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Table  1 


(Option  A)  cofferdam  removed  a.  Tended  to  produce  better  D.O.  distri¬ 
bution  during  early  phases  of 
stratification. 


b.  Resulted  in  most  favorable  D.O.  dis¬ 
tribution  in  late  periods  of  stratifi¬ 
cation  as  compared  with  other 
geometries. 

(Option  B)  cofferdam  left  a.  D.O.  distribution  tended  to  be  in¬ 

ferior  to  that  of  Option  A  during 
early  period  of  stratification. 

b.  Less  favorable  D.O.  distribution  in 
the  hypolimnion  during  the  late 
phases  of  stratification  (as  compared 
with  Option  A) . 

(Option  C)  breached  cofferdam  a.  Predicted  D.O.  distribution  identical 

to  unbreached  case. 

b.  Higher  potential  for  pulses  of  lower 
D.O.  releases  than  would  occur  with 
unbreached  cofferdam  (Option  B) . 

(Option  D)  raised  cofferdam  a.  Produced, anaerobic  conditions  in  the 

hypolimnion  between  June  and  mid- 
July.  Subsequently,  the  thickness  of 
the  anaerobic  zone  increased  to  be¬ 
tween  65  and  80  ft  by  the  end  of 
September. 


APPENDIX  A:  DISCUSSION  OF  NUMERICAL  MODEL 


Simulation  Model  Description 


1.  The  downstream  release  characteristics  and  internal  structure 
of  temperature  within  a  reservoir  were  predicted  with  a  numerical  simula¬ 
tion  model.  The  model,  hereinafter  identified  as  WESXEX*  was  devel¬ 
oped  at  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  based 
on  results  of  Clay  and  Fruh  (1970),  Edinger  and  Geyer  (1965),  Dake  and 
Harleman  (1966),  Bohan  and  Grace  (1973),  Imberger,**  and  others. 


Introduction 


2.  The  reservoir  is  conceptualized  as  a  number  of  homogeneous  hori¬ 
zontal  layers  stacked  vertically,  and  the  heat  sources  and  sinks  to  a  gen¬ 
eral  layer  are  represented  as  shown  in  Figure  Al.  The  solution  for  the 
temperature  history  of  a  general  layer  (or  control  volume)  is  obtained  by 
solving  a  conservation  of  mass  and  energy  equation.  The  governing 
equation  is: 


where 


30  Qi  9i  Qo  9o  1  _9  .  96  1  9(Qv9)  1  9H 

at  AAZ  "  AAZ  A  3Z  3Z  “  A  9Z  pC  A  3Z 

P 

Q  *  temperature  of  layer,  °F 
t  =  time,  days 

0.  =  inflow  temperature,  °F 

1  3 

Q.  *  flow  rate  into  layer,  ft  /day 

1  2 
A  =  horizontal  cross-sectional  area,  ft 

AZ  =  layer  thickness,  ft 

0  =  outflow  temperature,  °F 

°  3 

Qo  =*  outflow  rate,  ft  /day 

Z  =  elevation,  ft 


(Al) 


*  Bruce  Loctice,  "WESTEX — A  reservoir  Heat  Budget  Model”  (first  draft, 
unpublished),  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE, 
Vicksburg,  Miss. 

**  Private  communication,  Imberger,  and  Jorg. 
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k  -  vertical  diffueion  coefficient,  ft  /day 

3 

»  not  vertical  flow  into  or  out  of  layer,  ft  /day 
p  ■  density  of  water,  lb/ft^ 

Cp  ■  specific  heat  of  water,  Btu/lb/°F 
3H/3Z  -  external  heat  source,  Btu/ft/day 
Appropriate  boundary  conditions  must  be  supplied  at  the  air-water  and 
the  mud-water  interfaces.  Solution  of  Equation  1  for  each  layer  through 
time  yields  the  dynamic  vertical  temperature  distribution  of  the 
reservoir. 

3.  Conservative  water  quality  constituents  are  those  parameters 
for  which  the  concentration  in  an  elemental  volume  is  affected  only  by 
dilution,  not  by  biological  processes,  chemical  reactions,  or  physical 
phenomena  such  as  gravity.  Clay  and  Fruh  (1970)  demonstrated  that  suit¬ 
able  conservative  constituents  for  simulation  might  include  chloride, 
sulfate,  sodium,  or  silica  concentrations.  Conservative  constituent 
concentrations  can  be  budgeted  in  a  manner  analogous  to  that  used  for 
the  heat  budget.  Using  conservation  of  mass,  the  change  with  respect  to 
time  of  conservative  constituent  concentrations  can  be  represented  by: 


ic  «iCl  <!oCo  1  ,  1  8  >c\ 

St  “  A42  “  AAZ  ’  ~  A  3Z  A  K  I  II  I 

where 

C  *  concentration  of  a  conservative  constituent,  mg/i 
t  *  time,  days 

C.  *  inflow  concentration,  mg/fc 

1  3 

Q  *  inflow  rate,  ft  / day 

1  2 
A  »  cross-sectional  area,  ft 

&Z  ■  layer  thickness,  ft 

C  -  outflow  concentration,  mg /£ 

°  3 

Q  ■  outflow  rate,  ft  /day 

°  3 

-  net  vertical  flow,  ft  /day 
Z  ■  elevation,  ft 

2 

k  “  vertical  diffusion  coefficient,  ft  / day 
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4.  Bella  (1970),  Fruh  and  Davis  (1972),  and  Markofsky  and  Harle- 
man  (1971)  have  demonstrated  the  application  of  a  one-dimensional  model 
to  describe  a  nonconservative  constituent  budget.  Nonconservative  param¬ 
eters  can  be  budgeted  similar  to  conservative  parameters  if  the  noncon¬ 
servative  processes  are  appropriately  modeled.  Bella's  (1970)  develop¬ 
ment  of  a  one-dimensional  model  for  dissolved  oxygen  (D.O.)  structure  of 
a  reservoir  yielded  a  governing  equation  similar  to  Equation  2  with  an 
additional  terra  representing  the  result  of  photosynthetic  and  respira- 
tiou  processes  (source  and  sink  terms).  Fruh  and  Davis  (1972)  used  one- 
dimensional  models  to  describe  vertical  D.O.  structures  in  a  reservoir. 
The  initial  simulation  of  nonconservative  processes  with  WESTEX  was 
performed  by  Fontane  and  Bohan  (1974)  to  analyze  D.O.  in  Hartwell  and 
Clarks  Hill  Reservoirs.  Analysis  of  the  D.O.  distribution  was  accom¬ 
plished  in  this  study  in  a  manner  generally  ronsiatent  with  the  work  of 
Fontane  and  Bohan  (1974).  The  upper,  well-mixed  layers  in  the  reservoir 
were  assumed  to  be  approximately  saturated  (80  percent)  with  D.O.  This 
degree  of  saturation  wa3  assumed  to  extend  to  the  depth  at  which  a  tem¬ 
perature  difference  of  1°C  from  the  surface  temperature  existed.  A 
diffusion  analogy  similar  to  that  for  thermal  diffusion  (described  in 
paragraph  15)  was  used  to  obtain  the  desired  downward  transport  of  D.O. 
below  this  depth.  A  constant  daily  D.O.  depletion  rate  of  0.1  mg/i  was 
then  applied  to  each  layer  below  this  depth  in  the  reservoir.  The 
selection  of  the  1.0°C  temperature  difference  depth  and  the  0.1  mg/L 
daily  depletion  rate  was  based  on  an  analysis  of  observed  temperature 

and  D.O.  profiles  in  Hartwell  and  Clarks  Hill  Reservoirs  by  Foncane 
and  Bohan  (1974), 

5.  Fundamental  assumptions  and  the  various  processes  addressed 
in  the  mathematical  model  will  be  discussed  in  the  following  sections. 

Fundamental  Assumptions 

6.  Reservoir  hydrodynamic  phenomena  and  a  t\i,. renal  energy  balance 
are  used  to  predict  temperature  profiles  and  release  temperatures  in  the 
time  domain.  The  model  includes  computational  methods  for  simulating 
heat  transfer  at  the  air-water  interface,  advective  heat  due  to  inflow. 
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ouflow,  and  pumpback  processes,  and  the  internal  dispersion  of  thermal 
energy.  The  model  is  conceptually  one-dimensional  based  on  the  division 
of  the  impoundment  into  discrete  horizontal  layers  of  uniform  thickness. 
Assumptions  include  the  following: 

,a.  Isotherms  are  parallel  to  the  water  surface  both  laterally 
and  longitudinally. 

b.  The  water  in  each  discrete  layer  is  physically  homogeneous. 

c.  Internal  ndvection  and  heat  transfer  occur  only  in  the 
vertical  direction. 

d>  External  advection  (inflow,  outflow,  and  pumpback)  occurs 
as  a  uniform  distribution  within  each  layer. 

ji.  Internal  dispersion  (between  layers)  of  thermal  energy  is 
accomplished  by  a  diffusion  mechanism  that  combines  the 
effects  of  molecular  diffusion,  turbulent  diffusion,  and 
thermal  convection. 

7..  The  surface  heat  exchange,  internal  mixing,  and  advection  pro¬ 
cesses  are  simulated  separately  and  their  effects  are  introduced  sequen¬ 
tially  at  time  intervals.  A  simplified  flow  chart  of  the  mathematical 
simulation  procedure  is  presented  in  Figure  A2.  Each  process  will  be 
discussed. 

Surface  Heat  Exchange 

8.  The  net  heat  exchange  at  the  surface  is  composed  of  seven 
heat  exchange  processes: 

a. .  Shortwave  solar  radiation. 

b.  Reflected  shortwave  radiation. 

£.  Long-wave  atmospheric  radiation. 

d.  Reflected  long-wave  radiation. 

Heat  transfer  due  to  conduction. 

£.  Back  radiation  from  the  water  surface. 

£.  Heat  loss  due  to  evaporation. 

9.  The  surface  heat  transfer  process  is  solved  in  the  WESTEX 
model  with  an  approach  developed  by  Edinger  and  Geyer  (1965) .  The  ther¬ 
mal  equation  quantifying  the  net  surface  heat  exchange  (after  some 
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REPEAT  FOR  EACH  DAY 


Figure  A2.  Simplified  WESTEX  flow  chart 


linearization)  is: 


H  -  K  (E  -  es) 


(A3) 


where 

2 

H  ■  net  rate  of  surface  heat  transfer,  Btu/ft  /day 

2 

K  ■  coefficient  of  surface  heat  exchange,  Btu/ft  / day/°F 

E  ■  equilibrium  temperature,  °F 

9  ■  surface  temperature,  °F 

d 

Equilibrium  temperature  is  defined  as  that  temperature  at  which  the  net 
rate  of  heat  exchange  between  the  water  surface  and  the  atmosphere  is 
zero.  The  coefficient  of  surface  heat  exchange  is  the  rate  of  heat 
transfer  at  the  air-water  interface.  The  computation  of  equilibrium 
temperature  and  heat  exchange  coefficient  is  based  solely  on  meteorologi¬ 
cal  data  as  outlined  by  Edinger,  Duttweiler,  and  Geyer  (1968). 

10.  Heat  from  radiation,  with  the  exception  of  shortwave  radia¬ 
tion,  is  immediately  absorbed  at  the  surface  or  the  top  few  feet.  De¬ 
pending  upon  the  color  and  clarity  of  the  water,  shortwave  radiation 
penetrates  to  and  increases  the  temperature  at  greater  depths.  Based  on 
laboratory  investigations,  Dake  and  Harleman  (1966)  have  proposed  an 
exponential  decay  with  depth  for  describing  the  heat  flux  due  to  short¬ 
wave  penetration.  This  approach  is  used  in  WESTEX. 

11.  The  surface  heat  exchanges  are  implemented  in  the  model  by 
the  placement  of  a  percentage  of  the  incoming  shortwave  radiation 
throughout  the  depth  of  the  lake  and  by  the  placement  of  all  other 
sources  of  surface  heat  exchange  into  the  surface  layer.  The  shortwave 
radiation  is  distributed  exponentially,  so  that  most  is  absorbed  in  the 
top  layers  and  relatively  little  is  absorbed  in  the  lower  layers.  The 
procedure  can  be  expressed  mathematically  by  the  following  two  equations: 

H  -  K(E  -  0)  -  (1  -  B)x  (AA) 

s 

Hi  -  (1  -  8)  xe  AZi  (A5) 
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where 


H  ■  heat  transfer  rate  into  or  out  of  the  surface  layer, 

8  Btu/ft^/day 

$  *  percentage  of  the  shortwave  radiation  absorbed  in  the  surface 
layer 

2 

X  ■  total  incoming  shortwave  radiation  rate,  Btu/ft  /day 

2 

«  rate  of  heat  absorption  in  layer  i,  Btu/ft  /day 
e  ■  natural  logarithmic  base  (2.7183) 

X  »  heat  absorption  coefficient,  ft-'*' 
z^  «  depth  below  surface,  ft 

12.  Equations  A4  and  A5  are  applied  to  the  temperature  profile 
once  during  each  one-day  time-step.  The  net  heat  exchange  rate  into 
each  layer  is  computed  and  converted  into  a  temperature  change.  The 
temperature  changes  are  used  to  determine  an  updated  temperature  profile 
f o  r  the  lake  . 


Inf low 


13.  The  inflow  process  is  simulated  in  three  basic  steps.  The 
point  of  neutral  buoyancy  of  the  inflow  current  is  found,  water  in  the 
lake  is  displaced  and  mixed  upward,  and  the  water  In  the  lake  is  mixed 
with  the  inflow. 

14.  The  point  of  neutral  buoyancy  is  found  by  a  linear  interpola¬ 
tion  or  extrapolation  upon  the  density  profile  of  the  lake.  The  inflow 
volume  is  allocated  to  the  layer  of  neutral  buoyancy.  Where  inflow 
volume  into  a  layer  (i)  exceeds  the  physical  capacity  of  the  layer  (i), 
the  excess  volume  is  added  to  the  volume  of  the  next  higher  layer 

(i  +  1).  The  excess  volume  from  layer  i  (which  is  added  to  layer 
(i  +  1))  is  mixed  with  the  volume  in  layer  (i  +  1).  If  the  volume  in 
layer  (i  +  1)  after  mixing  is  greater  than  the  physical  capacity  of 
layer  (i  +  1),  the  excess  mixed  volume  is  placed  in  the  next  layer 
(i  +  2).  The  sequential  mixing  process  is  repeated  until  the  inflow  is 
distributed  within  the  water  column. 
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Internal  Mixing 


15.  The  internal  dispersion  process  is  represented  by  a  mixing 
scheme  mathematically  similar  to  a  simple  diffusion  analogy.  The  mixing 
process  is  used  to  account  for  the  effects  of  molecular  diffusion,  tur¬ 
bulence  diffusion,  and  additional  internal  mixing  processes  which  are  not 
explicitly  addressed.  Internal  mixing  transfers  heat  and  other  water 
quality  constituents  (e.g.,  D.O.)  between  adjacent  layers.  The  magni¬ 
tude  of  the  transfer  between  two  layers  in  the  WESTEX  model  is  expressed 
as  a  percentage  of  the  total  transfer  required  to  mix  the  two  layers. 

This  percentage  is  a  mixing  coefficient  that  is  defined  for  every  layer. 
Data  input  includes  values  of  the  mixing  coefficient  at  the  maximum  pool 
elevation  (a^)  and  at  the  bottom  of  the  lake  (c^) •  The  values  are  pred¬ 
icated  upon  field  data  if  available.  The  exponential  fit  between  the 
two  extreme  values  is  used  to  determine  the  appropriate  coefficient  for 
each  layer.  The  same  coefficients  are  used  for  mixing  qualities  as  for 
mixing  heat. 


Outflow 


16.  The  outflow  component  of  the  model  incorporates  the  selective 
withdrawal  techniques  for  orifice  flow  developed  at  WES  (Bohan  and  Grace 
1973).  Additionally,  the  selective  withdrawal  description  was  analyzed 
to  include  the  appropriate  application  to  simultaneous  releases  from 
multiple  horizontal  ports.  Transcendental  equations  define  the  location 
of  zero-velocity  limits  and  are  solved  with  a  half-interval  search 
method.  The  location  of  the  zero-velocity  limits  is  functionally  depen¬ 
dent  on  the  geometric  configuration  of  the  withdrawal  device,  release 
flow  rates,  and  density  structure  within  the  lake.  W^th  knowledge  of 
the  limits  of  the  withdrawal  zone,  the  velocity  profile  produced  by  the 
outflow  can  be  determined.  The  flow  from  each  layer  is  then  computed 
as  the  product  of  the  velocity  in  the  layer  and  the  width  and  thickness 
of  the  layer.  A  flow-weighted  average  is  applied  to  the  temperature  and 
D.O.  profiles  to  determine  the  value  of  the  release  temperature  and 
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D.O.  concentration  for  each  flow  condition. 


Operation  Schedules 

17.  Because  the  generation  and  pumpback  flow  rates  vary  consider¬ 
ably  during  a  day's  operation  and  these  flow  rates  affect  withdrawal  and 
pumpback  characteristics,  it  was  not  adequate  to  use  daily  average  flows 
to  accurately  simulate  generation  and  pumpback  operations.  Dally  aver¬ 
age  flow  routings  usually  are  quite  accurate  for  less  dynamic  reser¬ 
voirs  such  as  water  supply  and  flood-control  projects.  For  this  model, 
the  planned  generation  and  pumpback  flow  rates  and  durations  (in  hours) 
are  used  rather  than  daily  averaged  flow  rate3.  The  model  can  handle 

a  variety  of  different  flow  conditions  within  a  simulation  day.  Input 
includes  modes  of  operation  (generation  or  pumpback),  flow  rates,  and 
durations  of  flow  for  each  day. 

Pumpba* k  Characteristics 

18.  For  numerical  simulations  it  is  necessary  to  define  the  tem¬ 
perature  and  D.O.  content  of  the  flow  pumped  from  the  tailrace  of  Clarks 
Hill  Lake  into  RBRR.  The  pumpback  operation  was  modeled  by  mixing  the 
release  volume  from  RBRR  with  an  equivalent  surface  volume  from  Clarks 
Hill.  It  should  be  noted  that  pumpback  flows  containing  25  and  75  percent 
(by  volume)  of  Clarks  Hill  water  were  also  simulated.  In  this  range,  the 
fraction  of  pumpback  (assumed  to  be  contributed  from  Clarks  Hill)  had 
little  impact  on  the  release  water  quality.  The  temperature  and  D.O. 

of  the  contribution  from  RBRR  is  known  from  previous  generation  cycles. 

The  surface  temperature  for  Clarks  Hill  is  assumed  to  be  the  surface 
temperature  of  RBRR.  The  D.O.  content  of  the  contribution  from  Clarks 
Hill  is  saturation  corresponding  to  the  surface  temperatures  of  RBRR. 

19.  The  physical  model  (described  in  Appendix  B)  was  used  to 
determine  the  source,  magnitude,  and  distribution  of  the  reservoir 
waters  entrained  into  the  pumpback  jet.  Removal  of  the  entrained  volume 
from  the  reservoir  is  mathematically  simulated  in  a  manner  analogous  to 
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the  outflow  process,  except  that  the  outflow  quantity  and  quality  are 
mixed  with  the  pumpbaclt  volume.  Placement  and  mixing  In  the  reservoir 
of  the  augmented  pumpback  volume  are  analogous  to  the  Inflow  process 
described  in  paragraph  13. 

Density  Stability 

20.  The  WESTEX  model  requires  a  stable  density  profile  prior  to 
outflow.  Stability  is  assured  by  searching  adjacent  layers  from  bottom 
to  top  comparing  densities.  If  a  density  instability  is  identified,  the 
two  unstable  layers  are  partially  mixed,  and  the  mixed  density  is  com¬ 
pared  with  the  density  of  the  layer  above  the  mixed  region.  If  an  in¬ 
stability  still  remains,  the  layer  above  the  mixed  region  is  included  in 
the  mixed  region;  and  the  process  continues  until  stability  is  achieved 
or  the  surface  is  reached.  By  mixing  layers  above  an  instability  it  is 
possible  to  create  an  instability  below  the  mixed  region.  If  such  an 
instability  is  detected,  then  mixing  proceeds  downward  until  stability 
is  achieved  or  the  bottom  is  reached. 
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APPENDIX  B:  DISCUSSION  OF  NEAR-FIELD  PHYSICAL  MODEL  INVESTIGATION 


Purpose 

i 

i.  An  undistorted-scaLe  three-dimensional  near-field  model  of  the 
RBRR  project  was  constructed  in  order  to  evaluate  the  effect  of  va\ ious.: 
cofferdam  configurations  upon  the  withdrawal  and  pumpback  characteris- 
|  tics  of  the  project.  This  information  was  used  in  conjunction  with  a 

numerical  simulation  model  to  determine  which  design  option  was  most 

I 

i  beneficial  in  meeting  release  requirements  for  temperature  and  D.O.  at 

the  project. 

Model  Description 

i 

t 

i  2.  The  physical  model  (Figure  Bl)  was  constructed  to  a  1:150 

'  scale  (model  to  prototype)  and  reproduced  the  area  of  the  project  be- 

]  tween  the  dam  and  a  point  approximately  3000  ft  upstream  (shown  in 

i 

Figure  B2).  Prototype  penstock  intakes  (Figure  B3)  were  also  modeled 
;  upscream  of  the  transition  from  a  circular  to  a  rectangular  cross 

i  section.  Topography  was  reproduced  with  shaped  urethane  foam,  plywood, 

sheet  metal,  and  clear  acrylic  plastic.  The  cofferdam  was  constructed 
|  of  sheet  metal  to  simplify  modifications  to  the  model  for  each  cofferdam 

configuration.  Operation  of  each  model  penstock  was  controlled  by  a 
i  ball  valve.  These  valves  were  either  fully  opened  or  closed  during  with- 

!  drawa1  and  were  adjusted  to  equally  distribute  flow  through  the  four 

units  operating  during  pumpback.  Units  .1  through  4  are  reversible  tur- 
I  bines;  the  remaining  units  5  through  8  were  constructed  similarly  but 

I  without  the  capability  for  pumpback  operations. 

I 

j  Experimental  Procedure 

i 

3.  During  selective  withdrawal  testing  of  all  options,  the  out- 

%  _ 

?  flow  through  each  penstock  was  held  constant  at  7500  cfs.  Pumpback 

flow  was  constant  at  20,000  cfs  and  divided  equally  between  the  four 
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reversible  units.  Variations  in  total  outflow  rates  were  determined 
by  the  number  of  penstocks  withdrawing  from  the  reservoir.  The  major¬ 
ity  of  the  tests  were  conducted  with  two,  four,  or  six  penstocks  in 
operation. 

A.  Stratification  was  established  in  the  model  using  fresh  and 
saline  waters  to  simulate  the  density  distribution  expected  in  the  pool 
during  the  summer  stratification  period.  An  initial  density  profile 
was  determined  and  outflow  was  initiated  at  the  desired  test  flow. 
Vertical  dye  streaks  were  used  to  observe  the  formation  of  the  with¬ 
drawal  zone  and  to  determine  when  steady-state  conditions  were  estab¬ 
lished.  After  establishment  of  steady-state  conditions,  displacement 
of  the  vertical  dye  streaks  was  recorded  on  a  video  recorder  and  sub¬ 
sequently  the  velocity  distribution  and  limits  of  withdrawal  were 
determined. 

5.  Entrainment  was  determined  for  the  pumpback  tests  by  measuring 
the  dilution  of  the  pumpback  jet  with  a  fluorescent  dye.  The  epilimnion 
or  the  hypolimnion  was  dyed  in  separate  tests  and  the  contribution  of 
each  to  the  total  pumpback  volume  was  measured.  The  total  pumpback  vol¬ 
ume  consisted  of  the  volume  of  epilimnion  and  hypolimnion  waters  en¬ 
trained  plus  the  volume  of  the  pumpback  jet.  Additionally,  dye  streaks 
were  used  to  indicate  the  distribution  of  flow  within  the  model  during 
the  pumpback  operation. 


Testing 

6.  Model  tests  were  divided  into  four  options  with  each  option 
representing  a  different  upstream  cofferdam  configuration.  Option  A 
simulated  removal  of  the  upstream  cofferdam  prior  to  filling  the  lake. 
Option  B  simulated  the  cofferdam  as  initially  constructed  with  crest 
at  el  378.  In  Option  C,  the  Option  B  design  was  breached  to  simulate 
the  planned  second-stage  diversion.  In  Option  D,  the  same  cofferdam 
described  in  Option  B  was  modified  with  a  sheet-pile  wall  to  raise  the 
crest  to  el  400. 
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Analyst 8 


7.  In  the  analysis  of  model  data,  the  penstocks  were  treated  as 
orifices  and  the  cofferdam  was  treated  es  a  submerged  weir.  J.n  many 
cases,  the  analysis  was  conducted  using  the  same  dimensionless  group¬ 
ings  employed  by  Bohan  and  Grace.  They  experimentally  investigated  the 
stratified  flow  field  produced  by  withdrawal  through  orifices  (circular 
and  rectangular)  and  over  submerged  weirs.  The  weir  study  focused  on 
flow  over  a  weir  in  which  the  thermocline  was  below  the  crest  of  the 
weir.  Bohan  and  Grace  determined  dimensionless  groupings  that  char¬ 
acterized  the  flow  regime  and  subsequently  developed  semicmpirical 
descriptions  of  the  withdrawal  regime  for  the  condition  Investigated. 
These  descriptions  were  incorporated  into  the  numerical  routine  SELECT. 
However,  the  flow  conditions  in  RBRR  included  a  number  of  geometrical 
and  flow  conditions  for  which  SELECT  was  not  developed.  These  included: 

£.  Multiple  horizontal  ports. 

t>.  Bellmouthed  intake  geometries. 

c.  Submerged  weir  with  crest  elevation  below  the  elevation  of 
the  thermocline. 

Thun,  a  near-field  physical  model  was  constructed  to  acquire  the 
experimental  data  required  to  develop  a  description  of  the  unique  with¬ 
drawal  characteristics  applicable  tc  RBPft.  A  review  of  the  selective 
withdrawal  descriptions  developed  for  the  various  geometrical  configura¬ 
tions  will  be  discussed  in  four  sections.  They  are: 

a_.  Option  A:  cofferdam  removed. 

b.  Option  B:  cofferdam  left  in  place  (with  el  378). 

£.  Option  C:  cofferdam  left  in  place  with  a  second-p  ige 
diversion  breach. 

cL  Option  D:  cofferdam  crest  raised  from  el  378  to  el  400. 

8.  The  subsequent  discussion  of  the  results  obtained  from  Option  A 
includes  rn  overview  of  the  idealized  selective  withdrawal  work  of  Bohan 
and  Grace.  The  Bohan  and  Grace  descriptions  were  used  for  flow  condi¬ 
tions  which  the  SELECT  assumptions  were  valid.  Discussion  of  the 
specific  hydrodynamic  effects,  and  the  descriptions  characterizing  these 
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effects,  introduced  for  RBRR  will  be  given  thereafter.  Additionally,  at 
the  close  of  this  appendix  the  hydrodynamic  considerations  of  pumpback 
operation  are  discussed  for  each  option. 

Option  A 


i 

] 

i 


9.  In  Option  A,  the  cofferdam  was  removed  leaving  the  approach 
upstream  of  the  dam  unobstructed.  This  was  somewhat  similar  to  a  config¬ 
uration  studied  by  Bohan  and  Grace.  They  found  that  the  upper  and  lower 
limits  of  the  withdrawal  zone  should  be  described  as  a  function  of  the 
densimetric  Froude  number  such  that 


V 

o 


(Bl) 


where 

Vq  «  average  velocity  through  the  penstock,  fps 

Z  ■  vertical  difference  in  the  elevation  of  the  limit  in  question 
and  the  center  line  of  the  orifice,  ft 

2 

Aq  -  area  of  the  penstock,  ft 

Ap  -  density  difference  of  fluid  between  the  elevations  of  the 
center  line  of  the  orifice  and  the  limit  in  question,  g/cc 

p  -  density  at  the  elevation  of  the  center  line  of  the  orifice, 
g/cc 

2 

g  ■  acceleration  due  to  gravity,  ft/ sec 
The  quantity 


is  defined  as  the  densimetric  Froude  number,  where  gAp/p  represents 
the  gravitational  effect.  Bohan  and  Grace  found  that  this  function  was 
valid  for  both  upper  and  lower  limits  of  withdrawal,  except  in  cases 
where  the  upper  and  lower  limits  intersected  either  the  free  surface  or 
the  bottom  boundary,  respectively. 

10.  Bohan  and  Grace  predicted  the  location  of  maximum  velocity 
from  the  following  equation 
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where 

Yt  -  distance  from  lower  limit  to  location  of  maximum  velocity,  ft 

L 

H.  -  thickness  of  withdrawal  zone;  distance  from  lower  limit  to 
upper  limit,  ft 

Z  ■  distance  from  lower  limit  to  orifice  center  line,  ft 

They  found  that  the  velocity  profile  above  or  below  the  elevation  of 
maximum  velocity  could  be  calculated  for  boundary  interference  by: 


and  for  the  case  when  the  withdrawal  zone  is  not  limited  by  a  boundary: 
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where 

B  ■  2  for  orifice  and  weir  flow 
N  =  2  for  orifice  flow 
N  *  3  for  weir  flow 

V/V  x  *  rati°  of  local  velocity  to  maximum  velocity 

y  »  magnitude  of  distance  from  elevation  of  maximum  velocity 
to  some  local  elevation,  ft 

Y  ■  magnitude  of  distance  from  elevation  of  maximum  velocity 
to  the  withdrawal  limit  of  interest,  ft 


Ap  ■  density  difference  between  elevation  of  maximum  velocity 
and  local  elevation,  g/cc 


some  of  the  descriptions  in  SELECT  could  be  used  to  characterize  selec¬ 
tive  withdrawal  in  RBRR.  Equation  B1  was  used  to  predict  the  upper 
limits  and  Equation  B2  was  used  to  predict  the  point  of  maximum  velocity. 
A  new  description  for  the  lower  limit  of  withdrawal  was  required.  Analy¬ 
sis  of  the  observed  lower  limits  for  Option  A  yielded  the  following 
description: 


12.  Equations  B3  and  B4  were  applicable  for  Option  A.  Analysis 
of  the  data  indicated  that  the  appropriate  values  of  the  exponents  B 
and  N  were  1,40  and  4.15,  respectively. 

Option  B 

13.  Option  B  examined  the  hydrodynamic  effects  introduced  by 
leaving  the  cofferdam  in  place.  The  cofferdam  would  be  submerged  with 

a  nominal  crest  elevation  of  378.  The  elevation  of  the  thermocline  will 
be  above  the  weir  crest.  Previous  studies  have  primarily  addressed  with¬ 
drawal  characteristics  of  submerged  weirs  with  crests  above  the  eleva¬ 
tion  of  the  thermocline.  As  a  result  several  new  descriptions  had  to  be 
developed. 

14.  The  controlling  flow  structure,  either  the  weir  or  the  pen¬ 
stock,  was  found  to  be  different  for  each  of  the  two  limits.  The  func¬ 
tional  relationships  for  the  upper  limit  were  found  to  be  influenced 

by  penstock  effects;  those  of  the  lower  limits  were  affected  by  weir 
effects.  Thus,  no  single  functional  form  could  be  obtained  that  was 
adequate  for  the  description  of  both  limits.  The  equation  which  most 
adequately  described  the  position  of  the  upper  limit  of  withdrawal  for 
Option  B  was 


V  =  1.13 
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where 


H  ■  vertical  distance  from  bottom  of  lake  to  center  line  of 
0  orifice,  ft 

Z  -  vertical  distance  from  the  elevation  of  the  orifice  center 
line  to  the  upper  limit  of  withdrawal,  ft 

and  all  other  parameters  are  defined  in  Equation  Bl. 

15.  The  lower  limit  of  withdrawal  was  modeled  in  a  similar  man¬ 
ner  as  the  upper  limit.  The  densimetric  Froude  number  was  slightly  re¬ 
defined  for  weir  control  as 


where 

F  =  redefined  densimetric  Froude  number  for  weir  flow 
w 

Vw  *  average  velocity  over  the  weir,  fps 

Ap  *  density  difference  of  fluid  between  the  elevations  of  the 
w  weir  and  lower  limit,  g/cc 

P  =  density  of  fluid  at  the  elevation  of  the  weir,  g/cc 
w  2 

g  =  acceleration  due  to  gravity,  ft/sec 

■  vertical  distance  from  the  elevation  of  the  weir  crest  to  the 
lower  limit  of  the  withdrawal,  ft 

Analysis  of  the  Option  B  observed  lower  limits  resulted  in  the  following 
description 


where  Area  equals  area  of  flow  between  the  elevations  of  the  upper 

2 

limit  of  withdrawal  and  the  weir  crest,  ft  . 

16.  The  point  of  maximum  velocity  for  Option  B  was  defined  by 
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17.  Analysis  of  Option  B  test  results  indicated  that  Equations  B3 
and  BA  were  applicable  with  the  appropriate  values  of  exponents  B  and 
N  which  were  determined  to  be  2  and  1.275,  respectively.  With  these 
values.  Equations  B3  and  B4  were  used  to  describe  the  velocity  profile. 
Option  C 

18.  In  Option  C,  a  diversion  breach  with  the  bottom  at  el  310 
was  placed  on  the  east  side  of  the  cofferdam  with  a  crest  at  el  378. 
Comparison  of  the  outflow  densities  from  Options  B  and  C  indicated  that 
only  3  percent  of  the  total  outflow  was  withdrawn  through  the  breach. 
Consequently,  the  descriptive  equations  developed  for  Option  B  apply 

to  Option  C.  The  breach  was  numerically  modeled  as  an  orifice  with  a 
center  line  at  el  350  and  an  8500-sq-ft  area.  Except  for  the  breach, 
the  cofferdam  was  identical  to  that  simulated  in  Option  B. 

Option  D 

19.  In  Option  D,  the  crest  of  the  cofferdam  would  be  raised  from 
el  378  to  el  400  by  installing  sheet  piling.  Analysis  of  the  observed 
velocity  profiles  indicated  that  the  location  of  the  upper  limit  was 
controlled  by  flow  over  the  cofferdam  as  compared  with  control  by  the 
penstocks  with  the  Option  B  cofferdam  crest  at  el  378.  Analysis  of  the 
density  differences  and  distances  between  the  crest  and  the  upper  limit 
as  a  function  of  the  number  of  penstocks  operating  in  the  model  resulted 
in  the  equation 
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where 

Apwu  =  density  difference  of  fluid  between  the  elevations  of  the 
cofferdam  crest  and  the  upper  limit,  g/cc 

Z  m  vertical  distance  from  the  elevation  of  the  cofferdam  crest 

WU 

to  the  upper  limit,  ft 
Hw  *  height  of  the  cofferdam,  ft 
Np  =  number  of  penstocks  ir.  operation 
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Equation  B9  is  site  specific  and  was  developed  for  an  equal  to 

89  ft. 

20.  The  lower  limit  was  controlled  by  the  cofferdam  in  a  manner 
similar  to  that  seen  in  Option  B  tests.  As  a  result,  the  lower  limit 
was  described  by  Equation  B7. 

21.  The  location  of  the  point  of  maximum  velocity  was  a  function 
of  the  position  of  the  lower  limit.  The  relationship  developed  from  the 
data  was 
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This  equation  is  project  specific  and  may  not  be  applicable  to  other 
projects. 

22.  The  shape  of  the  boundary  unaffected  velocity  profile  deter¬ 
mined  for  Option  D  data  is  described  by  the  equation 
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23.  The  boundary  controlled  velocity  profile  description 
developed  for  Option  B  was  also  valid  for  Option  D. 

Pumpback  Results 

24.  The  proposed  operations  schedule  indicates  that  pumpback  will 
be  initiated  very  soon  after  the  end  of  the  afternoon  generation  period. 
Therefore,  withdrawal  currents  were  set  up  in  the  model  prior  to  simulat¬ 
ing  the  pumpback  operation.  Pumpback  flow  entered  the  model  through 
units  1  through  4  at  5000  cfs  per  unit.  The  density  of  the  pumpback  was 
equal  to  or  less  than  that  of  the  outflow,  depending  on  surface  water 
conditions  in  Clarks  Hill  Lake. 
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Option  A 

25.  Due  to  the  45-deg  inclination  of  the  penstocks  the  pumpback 
flow  was  directed  to  the  water  surface  and  subsequently  deflected  down¬ 
ward  in  the  lake  approximately  150  ft  in  front  of  the  penstocks.  After 
plunging,  the  pumpback  current  rose  tc  a  level  of  neutral  density  and 
spread  laterally  as  a  density  current  between  the  epilimnion  and 
hypolimnion  of  the  reservoir.  A  schematic  of  the  flow  pattern  is  pre¬ 
sented  in  Figure  B4.  Dilution  measurements  using  rhodamine-wt  fluores¬ 
cent  dye  indicated  that  the  final  volume  of  the  density  current  pro¬ 
duced  by  near-field  entrainment  was  2.5  times  the  pumpback  volume. 
Fluorescent  dye  measurements  and  dye  streak  displacement  observations 
indicated  that  96  percent  of  the  entrained  volume  came  from  below  the 
thermocline.  Although  the  velocity  did  not  vary  with  depth  in  the  en¬ 
trained  current  from  the  hypolimnion,  the  entrained  volume  varies  with 
depth  because  cross-sectional  area  of  flow  varies  with  depth. 

Option  B 

26.  The  pumpback  flow  was  similar  to  that  observed  in  Option  A 
tests  except  that  a  small  amount  of  the  pumpback  volume  was  laterally 
deflected  by  the  cofferdam  to  each  side.  Most  of  the  flow  traveled  up¬ 
stream  at  an  elevation  just  above  the  cofferdam  crest.  As  the  flow 
reached  a  point  approximately  500  ft  upstream  of  the  cofferdam,  the 
pumpback  jet  began  to  transform  to  a  density  current. 

27.  Most  of  the  entrained  flow  was  from  the  hypolimnion  as  in 
Option  A  with  the  entrained  flow  crossing  the  east  half  of  the  cofferdam. 
Additional  entrainment  from  the  hypolimnion  occurred  at  the  cofferdam 
where  the  pumpback  jet  entrained  flow  from  lower  layers.  Surface- 
entrained  flow  was  observed  to  come  from  west  of  the  cofferdam. 
Approximately  80  percent  of  the  entrained  volume  was  supplied  by  the 
hypolimnion  with  the  remainiig  20  percent  from  the  epilimnion.  Total 
entrainment  was  1.6  times  the  pumpback  volume. 

Option  C 

28.  Pumpback  flow  for  this  design  was  similar  to  Option  B  with 
one  exception.  Much  of  the  flow  entrained  from  the  hypolimnion  passed 
through  the  breach  in  the  east  end  of  the  cofferdam  which  formed  strong 


B13 


i.  OPTION  B  AND  C 


c.  OPTION  D 

Figure  B4.  A  schematic  of  the  pumpback.  jet  flow  pattern 
for  each  option  investigated 


lateral  currents  on  the  upstream  face  of  the  cofferdam.  As  a  result, 
little  or  no  entrainment  occurred  as  this  pumpback  jet  passed  over  the 
cofferdam.  Total  entrainment  for  this  design  is  1.5  times  the  pumpback 
volume  with  the  same  distribution  as  in  Option  B. 

Option  D 

29.  The  pumpback  flow  for  Option  D  was  significantly  different 
from  Options  A,  B,  and  C.  The  elevated  cofferdam  deflected  approximately 
40  percent  of  the  pumpback  flow  along  the  south  side  of  the  cofferdam 
and  inhibited  entrainment  from  the  hypolimnion.  In  addition,  the  pump¬ 
back  jet  was  significantly  weakened  prior  to  crossing  the  cofferdam,  re¬ 
sulting  in  a  pumpback  flow  that  crossed  the  cofferdam  as  a  density 
current  with  a  lower  boundary  at  el  400. 

30.  Entrainment  was  primarily  from  the  epilimnion  (85  percent) 
and  was  from  west  of  the  cofferdam.  Approximately  10  percent  of  the 
entrainment  was  from  the  hypo  1 '.inn ion  above  el  400  and  the  remaining 

5  percent,  from  below  el  400.  Total  entrainment  was  1.0  times  the  pump¬ 
back  volume. 
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